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Abstract 

The latest lattice results on kaon decays and mixing are reviewed. The discussion 
is focused on recent theoretical progress and new numerical calculations which appeared 
after Kaon 1999. 



1 Introduction 



In this talk a review of theoretical and numerical results on kaon decays, obtained with lat- 
tice QCD after Kaon 1999, is presented. Several issues of interest for the lattice approach 
can also be found in the talk by J. Donoghue at this Conference, to which I will refer when 
necessary. This contribution is being written just after the Lattice 2001 Conference []l]] 
where new intriguing results for e'/e by the CP-PACS §] and RBC § Collaborations have 
been presented. I added a discussion of these results to the material presented at Kaon 
2001. The main topics in this review are: i) the determination of the AI = 1/2 K — > tttt 
amplitudes; ii) the main contributions to e'/e due to the electropenguin operators (Qj and 
Qs in the standard notation) and to the strong penguin operator Qq. The most important 
novelties since Kaon 1999 are the following: 
Theoretical advances 

• Detailed and extended studies of the chiral expansion of several relevant amplitudes 
in finite and infinite volumes, in the quenched and unquenched cases, have been done 
in refs. 0-f§]- Further one-loop chiral calculations which are particularly important 
to control finite volume corrections and quenching effects in actual numerical simu- 
lations remain to be done. These calculations will also be very useful to extract the 
value of the physical K — ► tttt amplitudes from the lattice K — > tt matrix elements 
and to extrapolate both K —* tt and K — » tttt matrix elements, which are at present 
computed at unphysical values of the quark masses, to the physical point; 

• The importance of Final State Interaction (FSI) effects for AI = 1/2 transitions was 
first noticed in refs. || 1C]. More recently it has been emphasized by Bertolini, Eeg 



and Fabbrichesi [11| and Pallante and Pich jl^] that they may have large effects also 
for e'/e. Although the quantitative evaluation of FSI effects is still controversial [O ], 
there is a general consensus that, for a reliable calculation of kaon decay amplitudes, 
it is necessary to have a good theoretical control of FSI. This problem poses serious 
difficulties to lattice calculations based on the extraction of the K — ► tttt amplitudes 
from the K — ► it matrix elements using chiral perturbation theory (%PT), as done 
by the CP-PACS § and RBC Collaborations §; 



The Maiani- Testa no-go theorem [14| has prevented for a long time any attempt to 
directly calculate K — > tttt matrix elements in an Euclidean lattice. An important 
step towards the solution of this problem has recently been achieved by Lellouch and 
Liischer in ref. [jl5[ (denoted in the following by LL), who derived a relation between 
the K — > tttt matrix elements in a finite Euclidean volume and the physical kaon- 
decay amplitudes. Further investigation in this direction has been recently done in 
ref. [16] and exploratory numerical studies are just beginning. A discussion of these 



papers will be presented in the following. 
New numerical results 

• A lattice calculation of matrix element of the electro-magnetic operator (tt ^^!-^ ), 
which is relevant for the CP violating Kl — ► 7r°e + e _ decay, has been performed in 
ref. 0; 



New results for K°-K° mixing with Improved Wilson Fermions 1 1£ , ^] and Domain 



Wall Fermions (DWF) [p0| , 21] have appeared. With some differences, they confirm 
previous lattice results p^]. Calculations of (K°\Qi\K°) for all possible operators Qi 



which can mediate K°-K° mixing in extensions of the Standard Model have be done 



in refs. g3j. These results are very useful to put severe constraints on FCNC 



parameters of SUSY models [24|; 

The SPQcdR Collaboration has presented results for the amplitudes (-7r + 7r°|(54|if + ) 
and (tt7t\Q7 : s\K°) i=2, obtained by following the approach of and |l6| of com- 
puting directly the K — > irir amplitude || ^5|. The chiral behaviour of the matrix 
elements has also been studied in view of the extrapolation to the physical point and 
the necessary calculations of the relevant chiral loops are under way; 

In a first exploratory study, the SPQcdR Collaboration has also observed the first 
signal for {-kit\Q- \K°) I=0 and (7r7r|Q 6 |if ) /=0 p3j; 



• As mentioned before, CP-PACS and RBC have presented results for the CP conserv- 
ing AI =1/2 amplitude, Ao, an d for e'/e, obtained from K — ► n matrix elements 
using soft pion theorems. Both groups find a negative value, e'/e ~ — 5 x 1CT 4 . 
The possible origin of the discrepancy with the experimental results [^, 27] will be 
discussed in detail. 

In this short review a severe selection of arguments has been necessary. The reader, 
interested to have more information about lattice calculations of weak decays, including 
heavy flavours, may address either the nice reviews by L. Lellouch |22] and C. Bernard p^] 
at Lattice 2000, or the fortcoming Proceedings of the 2001 Lattice Conference [fflj. 



2 General framework 

Physical kaon weak decay amplitudes can be described with negligible error (of C(^ 2 /M^ / ), 
where fi is the renormalization scale) in terms of matrix elements of the effective weak 
Hamiltonian 

(ir7r\H AS=1 \K) , (1) 
written as combination of Wilson coefficients and renormalized local operators 

H AS=1 = -^^a(M)4(M). (2) 

» i 

The sum is over a complete set of operators, which depend on fi. In general there are 
12 four-fermion operators and two dimension-five operators: a chromomagnetic one and 
an electromagnetic one. In the Standard Model, the contribution of the dimension-five 
operators is usually neglected (SUSY effects may enhance the contribution of the chro- 
momagnetic operator p9[ ). The calculation of the matrix elements must be done non- 
perturbatively, and this is the role of the lattice, whereas the Wilson coefficients can be 
computed in perturbation theory. 

Wilson coefficients and renormalized operators 

For 7i^ s=1 the Wilson coefficients have been computed at the next-to-leading order in 



refs. |30|]-|3f|. The perturbative calculation is reliable provided that the scale // is large 
enough, fj, 3> Aqcd- In this respect calculations performed below the charm quark mass 
(m c ~ 1.3 GeV) are, in my opinion, suspicious. In fact, either /i <C m c , and then 



perturbation theory is questionable, or \i ~ m c , and then the effective three-flavour weak 
Hamiltonian (with propagating up, down and strange quarks) cannot be properly matched 
to the four-flavour theory (up, down, strange and charm) because of the presence of oper- 
ators of higher dimension which contribute at 0(/j, 2 /m 2 ). Another important remark is in 
order. Wilson coefficients and matrix elements of the operators Qi(fi) separately depend 
on the choice of the renormalization scale and scheme. This dependence cancel in physi- 
cal quantities, up to higher-order corrections in the perturbative expansion of the Wilson 
coefficients. For this crucial cancellation to take place, the non-perturbative method used 
to compute the hadronic matrix elements must allow a definition of the renormalized op- 
erators consistent with the scheme used in the calculation of the Wilson coefficients. The 
lattice approach satisfies this requirement. 

Matching of bare (divergent) lattice operators, Qi(a) to the renormalized ones is ob- 
tained by computing suitable renormalization constants Zik(fia) 

«4/=0,2(m) = (^\Qi{^)\K) 1=0,2 = z ik(va)(mr\Q k (a)\K) I=0:2 , (3) 

k 

where a is the lattice spacing. The ultra-violet (UV) problem, which deals with the con- 
struction of finite matrix elements of renormalized operators constructed from the bare 



lattice ones, has been addressed in a series of papers f3q|-|41[ and is, at least in principle, 
completely solved. The remaining difficulties are practical ones. On the one hand, the 
AI = 1/2 operators suffer from power divergences in the ultraviolet cutoff, 1/a. These 
divergences, that cannot be subtracted using perturbation theory can be eliminated 
by performing suitable numerical subtractions. The subtraction procedure, however, suf- 
fers from systematic uncertainties which are difficult to keep under control, see below. On 
the other hand, the perturbative calculation of the logarithmically divergent (or finite) 
Zik(iia) is rather inaccurate. Several non-perturbative methods have been developed in 
order to compute Zik(fJ.a) non-perturbatively [23], [^]-[46] and the uncertainties vary be- 



tween 1%, for the simplest bilinear operators, to IOh-25%, in the case of the four-fermion 
operators of interest. An accurate determination of Zik(fJ.a) for the AI = 1/2 operators 
is missing to date and more work is needed in this direction. 

Since numerical simulations are performed at finite values of the lattice spacing, a -1 ~ 
2-^4 GeV, another source of uncertainty in the determination of the matrix elements 
comes from discretization errors. They are of 0{A.qcdq), Q{\p\ a ) or 0(m q a), where 
p is a typical hadron momentum and m q the quark mass. The simplest strategy to 
reduce discretization effects consists in computing physical quantities at several values of 
the lattice spacing and then extrapolate to a — > 0. A different approach, pioneered by 



Symanzik and extensively studied on the lattice [46], is to reduce discretization errors 
from O(a) to 0(a 2 ) by improving the lattice action and operators. This method, which 
can also be combined with the extrapolation to a — ► 0, has been succesfully applied 
to the determination of the quark masses and of matrix elements of bilinear operators. 
A systematic study of the improvement of four fermion operators is still to be done. 



Note that with DWF [47] or overlapping Fermions [48] the errors are automatically of 
0(a 2 ) ||]. These formulations of lattice QCD are, however, much more demanding in 
terms of computing resources. Discretization errors correspond to the matching problems 
in effective theories with a low cutoff recently discussed by Cirigliano, Donoghue and 
Golowich In this respect the problem is softer for the lattice approach since numerical 
simulations are already performed at relatively large scales. 



K —> it and K — > nn matrix elements 



Two main roads have been suggested in the past in order to obtain A} =0 2 (m) : 

• Compute the K — > and K — > 7r matrix elements (0|Qj(/i)|if) and (7r|Qj(Ai)|^) 
and then derive (7T7t | <5i (a*) | -K") 7=0,2 using soft pion theorems |?6], |37]]. In this case 
the K — > -kit amplitudes can be evaluated only at the lowest order of the chiral 
expansion; 

• Compute directly (im\QiQi)\K) i=q^ U, EDI HI H- 



The main difficulty in the latter case is due to the relation between K — > tttt matrix 
elements computed in a finite Euclidean space-time volume and the corresponding physical 
amplitudes (the infrared problem). In the approach where the A l I=02 (/i) are extracted 

from the one-pion matrix elements, (ir\Qi(fj,)\K), the problem is the size of the corrections 
which relate matrix elements in the chiral limit to the corresponding physical amplitudes. 
These corrections are expected to be much larger for 1 = amplitudes, because of the 
larger FSI . Both the approaches have their advantages and drawbacks to which most 
of the following discussion will be devoted. 

The infrared problem arises from two sources: i) the unavoidable continuation of the 
theory to Euclidean space-time and ii) the use of a finite volume in numerical simulations. 
An important progress has been achieved by LL, who derived a relation between the lattice 
K — > tttt matrix elements in a finite volume and the physical kaon-decay amplitudes . 
An alternative discussion of boundary effects and the LL-formula, based on a study of the 



properties of correlators of local operators was presented in ref. [16|. In this approach the 
LL-formula is derived for all elastic states under the inelastic threshold, with exponential 
accuracy in the quantization volume. As a consequence the relation between finite- volume 
matrix elements and physical amplitudes, derived by Lellouch and Liischer for the lowest 
seven energy eigenstates, can be extended to all elastic states under the inelastic threshold. 
It can also been explicitly demonstrated how finite volume correlators converge to the 
corresponding ones in infinite volume. The derivation of JO]] is based on the property of 
correlators of local operators which can be expressed, with exponential accuracy, both as 
a sum or as an integral over intermediate states, when considering volumes larger than 
the interaction radius and Euclidean times < t ~ L. It can also been shown that it is 
possible to extract K — * tttt amplitudes also when the kaon mass, mx, does not match the 
two-pion energy, namely when the inserted operator carries a non-zero energy-momentum. 
Such amplitudes are very useful, for example, in the determination of the coefficients of 
operators appearing at higher orders in xPT, as illustrated by the numerical results for 
AI = 3/2 transitions presented in sec. ^. 

Let us sketch now the derivation of the result with an illustrative example which is 
not explicitly written in K6| . Consider the following Euclidean T-products (correlators): 

G(t,t K ) = (0\T[J(t)Qi(0)K^t K )]\0) , 

Git) = (0\T[J(t)J(0)]\0) , G(t K ) = (0\T[K(t K )K\0)\0) , (4) 

where J is a scalar operator which excites (annhilates) zero angular momentum tttt states 
from (to) the vacuum and K> is a pseudoscalar source which excites a kaon from the 
vacuum (t > and ij< < 0). At large time distances we have 

G(t,t K ) - V ^(0|J|7r7r(n)) y (vr7r(n)|Q J |iY)v(^|^ t |0)v eM~W n t - m K t K ) , 



G(t) -> V ^(0|J|7r^(n))y(7r^(n)|J|0)y exp(-W n i), 

-> F(0|^|K)y(K|^ t |0)y exp(-m^t^). (5) 
From a study of the time dependence of G(t, tx), G(t) and G(tx) we may extract 

1. the kaon mass m^', 

2. the two-pion energies on the finite lattice volume, W n . As shown by M. Liischer 
in J|l|], the W n are related to the infinite volume phase-shift of the two pions, S(k), 
via the following relations 

Wn =2f^, m±m =n , n=1 , 2 ,..., q = k ±, ( 6 ) 



where n is a non-negative integer [] and the function (j>{q) is defined in |51 

3. the operator matrix elements on a finite volume (inr(n)\Qi\K)y , (7T7r(n)| J\0)v and 
(K\K^\0) v . 

Moreover by a suitable choice of the lattice parameters it is possible to match one of the 
two-pion energies in such a way that mx = W n *. In practice, since it will be possible to 
disentangle only the first few states, the matrix elements will be computed with n* = 0-^2. 

The fundamental point is that it is possible to relate the finite- volume Euclidean matrix 
element {inv(n*)\Qi\K)y with the absolute value of the physical amplitude (Tm\Qi\K): 

(™\Qi\K) = ^{im(n*)\Qi\K) v , T = 32n 2 V 2 ^ E ^ K , ( 7 ) 

k[E) 




K(E) = \l—-ml, pv (E) = — = E. (8) 



where 



The last expression in eq. (|8[) is the one which one would heuristically derive from a naive 
interpretation of pv(E) as the density of states, cf. eq. (||). There are however, some 
technical subtleties with such an interpretation which will not be discussed here. Eq. (^) 
holds also when rax ^ W n and the operator carries non-zero energy-momentum |16| . 

By varying the kaon and pion masses and momenta, one may then fit the coefficients of 
the chiral expansion of (mr\Qi\K) and use these coefficients to extrapolate to the physical 
point. This procedure is necessary, at present, since it is not possible to compute the 
matrix elements at realistic values of the quark masses and in the unquenched case. To 
give an explicit example, the matrix elements of Q7 t s (f° r generic meson masses, one pion 
at rest and the other with a given momentum corresponding to an energy E n ) can be 
written as || 52] 



■ w , , i n ,'m K (m n + E n ) 
-iM 7 ^(m K , m^, E w ) = 7'' + d{ m w E n 



2 



+4' 8 iy- mK + — - m 7T E 7T j - 5l ,8 m K (m n + E w ) + (6 7 4 ' S + J^' 8 ) (2m 2 K + 4ml 
+5 7 6 ' 8 (4m 2 K + 2m 2 ^ + chiral logs + 0(p A ) (9) 



1 For n = 0, there are two solutions: one corresponding to k = which is spurious, the other giving the 
Liischer relation between the finite volume energy and the scattering length. 



A fit to the lattice data for M-j^, extracted from suitable correlation functions, allows the 
determination of the couplings 7 7,8 , (5^' 8 , etc. The results of the extrapolated amplitudes, 
if chiral logarithms are included, are independent of the cutoff used in the chiral theory. 
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Figure 1: Chiral behaviour of the matrix element of Q4, }A.^{mKi m -K- l E' V ). The extrapolation 
to the physical point using only the data corresponding to the lightest masses, without the effect 
of the chiral logarithms and using the operators of ref. Qj, is also shown as a shadowed band. 
The experimental result, computed as explained in the text, is also given as a line. 



3 Numerical Results 

In this section, the latest numerical results from lattice QCD for K — > tttt and K — ► 7r 
amplitudes and for Bk are reviewed, together with a comparison with other approaches. 

AI = 3/2 transitions and K°—K° mixing 

The SPQc*.dR Collaboration g, ^] has presented (quenched) results for both AI = 3/2 
and AI = 1/2 amplitudes, obtained from a calculation of K — > tttt matrix elements 
following the strategy of refs. |l5|, |lq| . For AI = 3/2 transitions, an extended study of 
the chiral behaviour of the matrix elements of the operators ^4,7,8, renormalized non- 
perturbatively, was performed. In figs, [j] and ||, M.i{mK,m 7T E 1T ) and M.^mK^m^E^) 
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Figure 2: Chiral behaviour of the matrix element of Qs, A4s{ m K, frim E^). The extrapolation 
to the physical point using only the data corresponding to the lightest masses, without the effect 
of the chiral logarithms and using the operators of ref. is also shown as a shadowed band. 



are shown as a function of the kaon mass. Note that the amplitude does depend on three 
independent quantities that are let to vary, namely mx, m n and E n (one of the two pions 
is always at rest). In the same figures the extrapolation to the physical point, performed 
by using the chiral expansion of the matrix elements as in eq. @, is given as a band. 
In fig. |l| the experimental number, extracted from the K + — > 7r + 7r decay rate using the 
Wilson coefficient of Q4 computed at the NLO, is also shown. The extrapolations are 
preliminary, since they do not include the effects of (quenched) chiral logarithms that 
have not been computed yet for the kinematical configurations used in this study. For 
some of the points, masses and momenta are too large to use chiral perturbation theory, 
and they have not been included in the fit. The preliminary results in fig. [l] already give 
us an interesting physics information. In the chiral limit, and using SU(3) symmetry, 
one may relate the K + — > 7r + 7r° amplitude to the K°-K° mixing parameter Bk [[5^]- 
In this limit, a large value for Bk, as found in lattice calculations and unitarity triangle 
analyses []54]], Bk ~ 0.85, would lead to a value of the K + — > 7r + 7r amplitude larger than 
the experimental one by ~ 50%. The extrapolation of the results of fig. [I], together with 
those for Bk obtained by the same collaboration [jl8| given in table ||, demonstrate that 
chiral corrections to both the K + — > 7r + 7r° amplitude and Bk can easily reconciliate a 
large value of Bk with the experimental K + — > 7r + 7r° amplitude. 

It is very instructive to compare the results for Q7 and Qs with those obtained in 
other approaches, table |]. For the sake of comparison all the results are converted to 
the MS scheme at a renormalization scale [i = 2 GeV. In the case of Qs we notice the 
nice agreement between lattice results obtained from K — > tt matrix elements using chiral 



perturbation theory and those from the first direct calculation of K — > tttt amplitudes. 
Similar numbers were obtained by the CP-PACS [§] and RBC Q collaborations, which 
computed K — > tt matrix elements with DWF. These groups presented the results in 
different renormalization schemes and scales and for this reason the values have not been 
included in the table. There is, however, a very large discrepancy of the lattice results 
with those obtained with dispersive methods [|55| or with the 1/N C expansion |5(|, whose 
results would correspond to a huge value of the S-parameter, B% = 3 -j- 7 R. In order to 
reproduce the experimental results for e'/e within the Standard Model, such a large value 
of £?8 implies a stratosferic value for B§. After Kaon 1999 a new analysis, performed with 
spectral function techniques by Bijnens, Gamiz and Prades appeared []58|j . In this paper, 
a value of (Qs) much lower than in refs. [55, 56] was found, although still about a factor of 
two larger than lattice determinations. The very low value of (Q7) found from the lattice 
K — > tttt calculation originates from large cancellations occuring in the renormalization 
of the relevant operator and requires further investigation. 



Table 1: K — > tttt matrix elements in GeV 3 , at fi = 2 GeV in the MS scheme, from lattice 
calculations (first three rows) and other approaches. 



Reference 


Method 




(Or) 


SPQcdR 
SPQcdR [ 

APE [§3 
Amherst [ 
Marseille [ 

BGP §8 


25] 2001 
lj 2001 
| 1999 
55 [ 1999 
56] 1998 
| 2001 


K — » TTTT 

K -> tt + xPT 
K -> tt + xPT 
Dispersion relations 

1/JVc 
Spectral functions 


0.53 ± 0.06 
0.49 ± 0.06 
0.50 ±0.10 
2.22 ±0.67 
3.50 ± 1.10 
1.2 ±0.5 


0.02 ±0.01 
0.10 ±0.03 
0.11 ±0.04 
0.16 ±0.10 
0.11 ±0.03 
0.26 ±0.03 



Lattice predictions for Bk have been very stable over the years, with a central value 
centered around 0.85. This is essentially also the value extracted from Unitarity Triangle 
Analyses [^4|], thus confirming that lattice QCD can predict (and not only postdict) 
physical quantities. In table |2| the lattice world average is given, together with some of the 
latest lattice results. Bk from CP-PACS |59| (with operators renormalized perturbatively) 
and RBC (2(J (with operators renormalized non-perturbatively) has been obtained with 
DWF. This formulation of QCD should garantee a better control of the chiral behaviour 
of the regularized theory with respect to Wilson-like fermions. The results from ref. [ |i"S| ] 
have been obtained using the non-perturbatively improved Wilson-like action, with a new 
method based on the Ward Identities |30| , which save us from the painful subtractions of 
the wrong chirality operators, which was necessary in the past. Similar strategies have 
been pursued with twisted-mass fermions [45|. The results with DWF are about 15% 



( 



below the world average and show a marked decrease at small quark masses. This could 
reconciliate the large value of Bk at the physical kaon mass with the low value of this 



parameter obtained in the chiral limit by ref. [31|. Whether the decrease at small quark 
masses is a physical effect, or is due to lattice artefacts (finite volume, residual chiral 
symmetry breaking etc.) remains to be investigated. 



2 Other results for Q-j and Qs can be found in refs. and |57 



Table 2: Bk in the MS scheme at the renormalization scale fi = 2 GeV and the renormalization 
group invariant Bk from recent lattice calculations. The results from SPQcdR, CP-PACS and 
RBC Collaborations are quenched. 



Reference 


Method 


5f s (2GeV) 


Br 


World Average 
SPQcdR |18| 2001 
SPQcdR M 2001 
CP-PACS gt 2001 

RBC [20| 2000 
Ciuchini et al. @] 2001 


with subtractions 
Ward identity method 

DWF Pert. Ren. 
DWF Nonpert. Ren. 
Unitarity Triangle 


0.63 ±0.03 ±0.10 
0.71 ±0.13 
0.70 ±0.10 
0.5746(61) (191) 
0.538 ± 0.08 
0.70l{u? 


0.86 ± 0.06 ± 14 
0.91 ±0.17 
0.90 ±0.13 
0.787 ± 0.008 
0.737 ±0.011 
Q0+ - 30 




Figure 3: The disconnected penguin (DP) diagram present in AI = transitions. 

For AI = 1/2 transitions and e'/e, the subtractions of the power divergencies, neces- 
sary to obtain finite matrix elements, are the major obstacle in lattice calculations. These 
divergencies arise from penguin contractions of the relevant operators, fig. ||, which induce 
a mixing with operators of lower dimensions, namely sa^ v G^t^d, s(j^ u ^^G^ v t A d, sd and 
s^d [37]. Power divergencies and subtractions are encountered with both the strategies 



(K — ► 7T7T or K — > ir) adopted to compute the physical amplitudes. For example, with a 



propagating charm quark, ReAo is computed in terms of the matrix elements of Q {Q\ 
and Q2) only. In this case, due to the GIM mechanism, the subtraction is implicit in the 
difference of penguin diagrams with a charm and an up quark propagating in the loop: 

Q ± = s7m (1 - 75) uu^ (1 - 75) d± S7 M (1 - 75) du^ (1 - 75) u - {u -»• c) . (10) 

In the past any attempt to compute AI = 1/2 K — > tttt matrix elements failed 
because no visible signal was observed after the subtraction of the power divergencies ]|6^, 



63]. These calculations, performed about twelve years ago with very modest computer 
resources, on small lattices and with little statistics, were abandoned after the publication 
of the Maiani- Testa no-go theorem |]14||. 
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Figure 4: First signal for the matrix element (7T7t|Q \K) at the matching point itlk ~ Wq, 
obtained with 400 configurations and V = 24 3 x 64. 



This year, the SPQc^R Collaboration has found the first signals for both the matrix 
elements of Q~ and Qq, figs. [| and ||. The difference with respect to previous attempts is 
given by a much larger volume and statistics, about 400 configurations with V = 24 3 x 64, 
compared to a few tens of configurations on lattices with V = 16 3 x 32 (8 configurations on 
V = 24 3 x 40) || or 110 configuration of a lattice with V = 16 x 12 x 10 x 32 (sic !!) pR . 
A crucial ingredient is also to work at the point corresponding to the matching condition 
nix = Wq, where Wq is the energy of the two pions at rest on the finite volume used 
in the simulation, namely with non degenerate strange and light quark masses []. Past 
calculations were always performed, instead, at the degenerate point, m# = m n . The 
results are very preliminary and in particular the two-pion energy in the finite volume for 
this channel gives a scattering lenght, a/=o, in disagreement with expectations both in 
value and in the mass dependence. This is to be contrasted with the AI = 3/2 case, where 
the analysis of the scattering length aj = 2 is in good agreement with expectations [f25|]. 
Much more work is needed to clarify these problems before trying to extrapolate the 
amplitude to the physical point. 

At the matching point mx = Wq, there are also results for (tttt\Qq\K) . In this case 
there is no GIM mechanism at work, and a finite subtraction of the matrix element of the 
pseudoscalar density operator Qp = (m s — m^s^d must be done (for parity, instead, the 
subtraction of the scalar operator Qs = (jn s + m^)sd is not necessary). The coefficient of 



In this case, as in the case of degenerate quark masses, the subtraction of operators of lower dimensions is 
not necessary [p9[. 
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Figure 5: Fzrsi signal for the matrix element (im\Q e \K) at the matching point uik ~ W , 
obtained with 400 configurations and V = 24 3 x 64. 
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Figure 6: (7T7r|s75(i|i^) at the matching point itlk ~ Wq, obtained with 400 configurations and 
V = 24 3 x 64. At the matching point (7T7r|s75(i|.K") decreases by a huge factor and thus the 
subtraction for (7T7t\Qq\K) is very small. 



the mixing of Qp with Qq is quadratically divergent f] 



Qe-Qe + ^^Qp (11) 

If not for the explicit chiral symmetry breaking of the lattice action, (irir\Qp\K) would 
vanish by the equation of motions when tuk = W n *. With an improved action, (irir\Qp\K) 
is of 0(a 2 ) and thus a finite subtraction is necessary. The preliminary results shown in 
fig. H show that that the value of (7T7r| s^ysdl-fC) drops by a factor of about 20 with respect to 
(7r7rlQ6l.fr) when rrix = Wo (it is of the same order when the kinematics is not matched) 
and thus the subtraction is not very critical. The signal itself is very noisy however, see 
fig. H|, and even with 400 configurations the statistical error on (tttt\Qq\K) is about 50%. 
Thus a rather large statistics will be necessary to obtain a relatively accurate result. 

K — > 7T7r amplitudes in the chiral limit can also be extracted from the calculation of 
K — > ir matrix elements, and indeed this has been the most popular method in lattice 
calculations of AI = 3/2 transitions. The main advantage of this approach is that the 
three-point correlators, necessary for the extraction of the matrix elements, are less noisy 
than the four-point correlators used in the K — > tttt case. Moreover, for K — > 7r matrix 
elements, finite volume corrections are exponentially small as L — > 00. The main disad- 
vantage is that the K — > tttt amplitudes are obtained, using soft pion theorems, in the 
chiral limit only and then the effect of FSI is definitively lost. If these are important for 
AI = 1/2 channels there is no hope, then, to recover the physical amplitudes. 

In order to make the K —* 7r matrix element finite at least a subtraction is nec- 
essary, even in the absence of explicit chiral symmetry breaking in the action. With 
Wilson-like Fermions, for which chiral symmetry is explicitly broken, the number of sub- 
tractions makes this approach unpracticable as also demonstrated by the failure of past 



attempts j62|, |63J. The method has acquired a renewed popularity, instead, with recent 
formulations of the lattice theory, DWF or overlapping Fermions, for which chiral sym- 
metry breaking is absent or strongly reduced in practice. In this case, and under the 
hypothesis that chiral symmetry is under control, a finite amplitude can be obtained by 
subtracting the scalar density amplitude with a suitable coefficient Cj which depends on 
the operator at hand 

(7r|Qf" b |^> = (w\Qi\K) - CMQs\K) . (12) 

The power divergent coefficients cannot be computed perturbatively. Following ref. [[3(| . 
when chiral symmetry is preserved, the Cj can be determined by the condition that the 
K — > matrix element of the subtracted operator vanishes 

(0\Qi -QQpIK) =0. (13) 

The coefficients Cj have been obtained either using non degenerate quarks, m s 7^ m^, 
by the RBC Collaboration, or from the derivatives of the 2-point correlation function 
with respect to the quark mass, by the CP-PACS Collaboration. Numerically, the matrix 
element of the subtracted operator is much smaller than the unsubtracted one, as shown 
by fig. [7], taken from the recent work of the RBC collaboration ||. This implies that any 



The mixing with the chromomagnetic operator is not discussed here. This mixing is a small, finite correction 
which can be ignored to simplify the discussion. 




Figure 7: (ir\Q\K) for the unsubtracted operator, Qq (squares), the subtraction, CqQs (Circles), 
and the total, Ql ub see eq. (f/l[), as a function of the quark mass irtf. The results are from the 
RBC Collaboration. 



systematic uncertainty which enters in the subtraction procedure can have huge effects in 
the determination of the physical amplitudes. After the subtraction, the chiral dependence 
of the matrix element has to be fitted in order to extrapolate it to the chiral limit. Both 
groups have included the logarithmic corrections which arise in quenched %PT in the fit 
and, in some cases, polynomial corrections of higher order in m\. The chiral behaviour 
observed by CP-PACS is very satisfactory, as shown by fig. ||, less good in the case of the 
RBC Collaboration. The difference may be due to the fact that a different gluon action 
is used in the two cases, corresponding to smaller chirally breaking effects for CP-PACS. 
This can be monitored by measuring the so called "residual mass" which should be zero 
for perfect chirality and is about a factor of ten smaller for CP-PCAS than for RBC. For a 
more extended discussion on this important point the reader can refer to the forthcoming 
proceedings of the Lattice 2001 Conference Q, or those of the 2000 edition 64]. 

Let me also mention that the overall renormalization constants have been computed 
perturbatively by CP-PACS [^0| and non-perturbatively by RBC ||44|| . My compendium 
of the physics results obtained by the two groups is given in table y 

A few observations are necessary at this point. First of all, ReAo, and consequently 
ReAo jReAi , from the RBC Collaboration are in good agreement with experimental val- 
ues, contrary to CP-PACS which finds ReAo smaller by about a factor of two than the 
experimental number. Since the two groups use the same lattice formulation of the the- 
ory and differ only by, hopefully, marginal details, the reason of this difference should 
be clarified. In both cases, however, e'/e is in total disagreement with the data. The 




Figure 8: Chiral behaviour of the subtracted matrix element (7i\Q^ ub \K) from the CP-PACS 
Collaboration. The curves represent different fits used to extrapolate to the chiral limit. 



Table 3: Lattice results for AI = 1/2 transitions using K —>■ n matrix elements from RBC and 
CP-PACS. The experimental numbers are also given. 



Reference 


ReAo 


ReA 2 


ReAo/ Re A 2 


e'/e 


RBC {§ 2001 
CP-PACS § 2001 
Exps. m, 27, H 2001 


29 -=- 31 x 10~ 8 
16 -=-21 x 10~ 8 
33.3 x 10" 8 


1.1 -=- 1.2 x 10~ 8 
1.3 1.5 x 10~ 8 
1.5 x 10~ 8 


24^-27 
9^ 12 
22.2 


-8-r-4x 10~ 4 
-7-=- -2 x 10~ 4 
17.2 ± 1.8 x 10" 4 



main reason is that the value of the matrix element of Qq is much smaller than what 
would be necessary to reproduce the experimental value (it corresponds approximatively 
to Bq = 0.3 ~r~ 0.4), see also the talk by J. Donoghue at this Conference. Let me list a 
number of sources of systematic errors which may explain these embarassing results: 

1. Chirality By working with DWF at a finite fifth dimension, iVg = 16, a residual 
chiral symmetry breaking is present in the theory. The amount of residual symmetry 
breaking is parametrized by a mass scale denoted by m res ~ 0.2 4- 2.0 MeV. In the 
presence of explicit chiral symmetry breaking, the coefficients Cj determined from 
K — ► differs from the correct one and this may induce an error of 0{m res a~ 2 ) ~ 
(200 MeV) 3 on matrix elements which are of the order of A-q CD ~ (300 MeV) 3 . Both 
groups claim to have this point fully under control 0, H], A calculation at a larger 
value of 2V5, with all the other parameters unchanged, would be very useful to clarify 
the situation. 

2. Matching below the charm mass Both groups have so far presented results at a 
renormalization scale just below the charm quark mass. As discussed in the in- 
troduction, the matching of the effective theory is rather problematic at such low 
scales. 



3. Extra Quenched Chiral Logarithms As shown by fig. [7|, the subtraction of the power 
divergencies is very critical. Besides the effects discussed in 1., there is another 
delicate problem which has been recently rised by Golterman and Pallante ||. In 
the quenched theory the (8, 1) operators, such as Q~ and Qq, do not belong anymore 
to irreducible representations of the chiral group and this gives rise to spurious chiral 
logarithms. These logarithms affect the subtraction procedure and have not been 
taken into account in the analyses of RBC and CP-PACS. 

4. Higher order chiral corrections and FSI Higher order chiral contributions, among 
which FSI have also to be accounted, can produce large variations of the matrix 
elements between the chiral limit and the physical point. In simulations where only 
K — > tt matrix elements are computed, K — > tttt amplitudes can only be obtained at 
lowest order in xPT where these physical effects are missing. 

5. New Physics As noticed by Murayama and Masiero [pSfl , it is possible to produce 
large effects for e'/e in SUSY without violating other bounds coming from FCNC. 
If the lattice results are correct, this is an open possibility. Before invoking new 
physics, at least the result for ReAo should be established with more confidence. 
Given that the two groups still don't agree on this quantity, more work is needed. 



4 Conclusion and outlook 

A renewed activity in lattice calculations of kaon decays and mixing has developed in 
the last two years. For AI = 3/2 transitions, K + — > 7r + 7r° and electroweak penguin 
amplitudes, new and more precise results have been obtained. In this case, removal of lat- 
tice artefacts, by extrapolation to the continuum and/or improvement, and unquenched 
calculations are around the corner. For AI = 1/2 transitions, it has been shown that 



direct computations of K —* 7T7t amplitudes, including FSI, is possible pq , 16], although 
the pratical implementation with sufficient accuracy will require more time. First results, 
obtained by using K — > tt lattice matrix elements computed with DWF and soft-pion 
theorems, show a striking disagreement with the experiments for e'/e. For ReAo, there 
is a factor of two between the values found by the CP-PACS @ and RBC || Collabora- 
tions. More work is needed in this case to clarify all these points. Many related physical 
quantities, like the strong interaction phase-shifts for tttt scattering, the chromomagnetic 
operator matrix elements, semileptonic K — > irirlvi amplitudes and the scalar semileptonic 
form factor will also be obtained as a byproduct from lattice investigations of non-leptonic 
kaon decays, see also J. Donoghue at this Conference. 
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